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ABSTRACT 

Three types of quasi-periodic oscillations (QPOs) have been discovered so far in 
the persistent emission of the most luminous neutron-star low-mass X-ray binaries, the 
Z sources: ~10-60 Hz horizontal and ~6-20 Hz normal/flaring branch oscillations, 
and ~200-1200 Hz kilohertz QPOs, which usually occur in a pair. Here we study the 
horizontal branch oscillations and the two simultaneous kilohertz QPOs, which were 
discovered using the Rossi X-ray Timing Explorer, comparing their properties in five Z 
sources with the predictions of the magnetospheric beat-frequency and Lense-Thirring 
precession models. 

We find that the variation of the horizontal-branch oscillation frequency with 
accretion rate predicted by the magnetospheric beat-frequency model for a purely 
dipolar stellar magnetic field and a radiation-pressure-dominated inner accretion disk 
is consistent with the observed variation. The model predicts a universal relation 
between the horizontal branch oscillation, stellar spin, and upper kilohertz QPO 
frequencies that agrees with the data on five Z sources. The model implies that the 
neutron stars in the Z sources are near magnetic spin equilibrium, that their magnetic 
field strengths are ~ 10^-10^° G, and that the critical fastness parameter for these 
sources is ^ 0.8. If the frequency of the upper kilohertz QPO is an orbital frequency 
in the accretion disk, the magnetospheric beat-frequency model requires that a small 
fraction of the gas in the disk does not couple strongly to the stellar magnetic field 
at 3-4 stellar radii but instead drifts slowly inward in nearly circular orbits until it is 
within a few kilometers of the neutron star surface. 

The Lense-Thirring precession model is consistent with the observed magnitudes 
of the horizontal branch oscillation frequencies only if the moments of inertia of 
the neutron stars in the Z sources are ~ 4-5 times larger than the largest values 
predicted by realistic neutron-star equations of state. If instead the moments of 
inertia of neutron stars have the size expected and their spin frequencies in the Z 
sources are approximately equal to the frequency separation of the kilohertz QPOs, 
Lense-Thirring precession can account for the magnitudes of the horizontal branch 
oscillation frequencies only if the fundamental frequency of the horizontal branch 
oscillation is at least four times the precession frequency. The trend of the correlation 
between the horizontal branch oscillation frequency and the frequency of the upper 
kilohertz QPO is not consistent with this model. 

We argue that the change in the slope of the correlation between the frequency of 
the horizontal branch oscillation and the frequency of the upper kilohertz QPO, when 
the latter is greater than 850 Hz, is directly related to the varying frequency separation 
of the kilohertz QPOs. 

Subject headings: accretion, accretion disks — stars: neutron — stars: individual 
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(Sco X-1, Cyg X-2, GX 5-1, GX 340+0, GX 17+2) — X-rays: stars 

1. INTRODUCTION 

The discovery in the Z-type low- mass X-ray binaries (LMXBs) of ~ 10-60 Hz quasi-periodic 
oscillations (QPOs) with centroid frequencies that arc positively correlated with mass accretion 
rate (van der Klis et al. 1985; see van der Klis 1989 for a review) has led to a significant 
improvement in our understanding of such systems. These horizontal-branch oscillations (HBOs), 
which are named after the branch in X-ray color-color diagrams where they appear, were the 
first rapid-variability phenomena discovered in LMXBs and have played a key role in organizing 
the complex phenomenology of these sources that has emerged over the past decade (see, e.g., 
Hasinger k, van der Klis 1989; van der Klis 1989). 

Very soon after the discovery of the HBO, its centroid frequency was identified with the 
difference between the Keplerian frequency at the radius where the neutron star magnetosphere 
couples strongly to the gas in the accretion disk and the spin frequency of the neutron star (Alpar 
&; Shaham 1985; Lamb et al. 1985; Shibazaki & Lamb 1987). This magnetospheric beat-frequency 
interpretation of the HBO was found to agree well with the observed properties of the HBOs, 
including the dependence of the HBO frequency on X-ray countrate (Alpar k, Shaham 1985; Lamb 
et al. 1985; Ghosh &; Lamb 1992), the existence of correlated low-frequency noise (Lamb et al. 
1985; Shibazaki &: Lamb 1987), and the absence of any detectable QPO with a frequency equal 
to the Keplerian frequency at the magnetic coupling radius (Lamb 1988). The magnetospheric 
beat-frequency model predicted that the neutron stars in the Z sources have spin frequencies 
~ 200-350 Hz and magnetic field strengths ~ 10^-10^° G (see Alpar & Shaham 1985; Ghosh 
& Lamb 1992; Wijnands et al. 1996), consistent with the hypothesis that these stars are the 
progenitors of the millisecond rotation-powered pulsars (Alpar &; Shaham 1985; see Alpar et 
al. 1982; Radhakrishnan & Shrinivasan 1982). The neutron star properties inferred from the 
magnetospheric beat-frequency model have subsequently been shown to be consistent with the 
magnetic field strengths inferred from models of the X-ray spectra of the Z sources (Psaltis, Lamb, 
k, Miller 1995; Psaltis k. Lamb 1998) and with the 290-325 Hz spin frequencies inferred from 
the frequency separation of the two simultaneous QPOs with frequencies ~ 1 kHz (hereafter the 
kilohcrtz QPOs) and from the high-frequency oscillations observed during type I X-ray bursts 
in other neutron-star LMXBs (Strohmayer et al. 1996; Miller, Lamb, k Psaltis 1998; Miller 
1999). If the upper kilohcrtz QPO is an orbital frequency in the inner disk, the magnetospheric 
beat-frequency model of the HBO requires that a small fraction of the gas in the accretion disk 
must penetrate to radii smaller than the radius where it initially couples to the stellar magnetic 
field (van der Klis et al. 1997; see Miller et al. 1998 for a discussion), because observations show 
that the kilohcrtz QPOs are present at the same time as the HBO (see, e.g., Wijnands k van der 
Klis 1998). In addition to the HBOs, the magnetospheric beat- frequency model has been used to 
explain successfully similar QPOs observed in accretion-powered pulsars, where the neutron star 
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spin frequency can be measured directly and the magnetic field strength can be estimated from 
the accretion torque, providing a stringent test of the model (see, e.g., Angelini, Stella, &; Parmar 
1989; Ghosh 1996; Finger, Wilson, & Harmon 1996). 

Recently, Stella & Vietri (1998) have proposed an alternative HBO mechanism, motivated 
by concern about whether orbiting gas can penetrate inside the magnetic coupling radius in the 
Z sources. In this model, the magnetic field of the neutron star plays no role in generating the 
HBO. Instead, the HBO observed in the Z sources and the power-spectral peaks with frequencies 
^ 20-60 Hz seen in some atoll sources (see, e.g.. Ford & van der Klis 1998) are both generated by 
nodal (Lense Thirring and classical) precession of a tilted ring of gas at a special radius in the 
inner disk. Stella &; Vietri suggested that the nodal precession frequency of the ring is visible in 
X-rays because of the changes in the Doppler shift of radiation from blobs orbiting in the ring, 
changes in occultations by such blobs, or the changing aspect of the ring seen by an observer. 
Subsequently, Markovic k. Lamb (1998) studied the normal modes of the inner disk and showed 
that typically ~ 10 high-frequency nodal precession modes are weakly damped. Nodal precession 
has also been proposed by Cui, Zhang, & Chen (1998; see also Ipser 1996) as an explanation for 
the QPOs observed in black hole candidates. If the HBO is generated by nodal precession at the 
same radius in the accretion disk where orbital motion generates the kilohertz QPO, as proposed 
by Stella & Vietri (1998), then the HBO frequency, the neutron star spin frequency, and the 
frequency of the upper kilohertz QPO should satisfy a specific relation. The shape of this relation 
was shown to be consistent with observations of the HBO and kilohertz QPO frequencies observed 
in the Z sources GX 17+2 and GX 5-1 (Stella 1997; Stella & Vietri 1997, 1998; Morsink & Stella 
1999), although the predicted precession frequencies were found to be smaller than the observed 
HBO frequencies. 

In this paper we use data on five Z sources obtained using the Rossi X-ray Timing Explorer 
(RXTE) to investigate further the origin of the HBO. All of these data have been fully reported 
elsewhere. The sources we consider are GX 17+2 (Wijnands et al. 1997; Homan et al. 1998), 
GX 5-1 (Wijnands et al. 1998b), GX 340+0 (Jonker et al. 1998), Cyg X-2 (Wijnands et al. 1998a), 
and Sco X-1 (van der Klis et al. 1997). In all of these sources the HBO and two kilohertz QPOs 
have been observed simultaneously. In GX 349+2, the sixth originally identified Z source, no 
HBO has so far been detected simultaneously with the kilohertz QPOs (Kuulkers & van der Klis 
1998; Zhang, Strohmayer, & Swank 1998). Therefore, we cannot include this source in the present 
study. We investigate the magnetospheric beat-frequency model in §2 and the Lense-Thirring 
precession model in §3, comparing their predictions with the available data. In §4 we summarize 
our conclusions and their implications for the properties of the neutron stars in the Z sources. 
Finally, we characterize the correlations between the various frequencies in a model-independent 
way in Appendix A, in order to facilitate comparison of the present data with future data or other 
theoretical models. 
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2. THE MAGNETOSPHERIC BEAT-FREQUENCY INTERPRETATION 

2.1. Model Predictions and Comparison with Observations 

In the magnetospheric beat-frequency model of the HBO (Alpar & Shaham 1985; Lamb et al. 
1985), the centroid frequency z^hbo of the HBO is identified with the beat between the Keplerian 
frequency z^k,™. at the radius where the neutron star magnetic field couples strongly to the gas 
in the accretion disk and the spin frequency Ug of the neutron star. The frequency of this beat is 

I'MBF = l^K,m. — ■ (1) 

Ghosh & Lamb (1992) computed the dependence of z^mbf on the stellar mass and magnetic 
moment and the accretion rate for a variety of simple models of the inner accretion disk. They 
found that if the coupling radius is in an asymptotic region of the disk, then 

I/MBF ^ VKflM'^^lj ( ^ ) -Vs , (2) 

where /X27 is the magnetic moment of the neutron star in units of 10^'' G cm'^, M is its gravitational 
mass in units of solar masses, M is the mass accretion rate, and Me is the Eddington critical 
mass accretion rate onto a neutron star of 10 km radius; the proportionality constant t'K,o and the 
exponents a, (3, and 7 in equation are different for different models of the inner accretion disk 
and are listed in Table 1. The dimensionless parameter ^ describes the fraction of the mass flux 
through the inner disk that couples to the stellar majgnetic field a,t t^^i and is introduced here to 
allow for the possibility that some of the gas in the disk does not couple to the stellar magnetic 
field at Vm but instead penetrates to smaller radii, as required if the upper kilohertz QPO is an 
orbital frequency in the inner disk (Miller et al. 1998); ^ may depend on the mass accretion rate. 

In the magnetospheric beat-frequency model of the HBO, the steep dependence of i^hbo on the 
mass accretion rate inferred from the EXOSAT data implies that i^mbf is small compared to VK,m 
and hence that the spin frequencies Vg of the neutron stars in the Z sources are very close to, but 
less than, VK,m (Alpar & Shaham 1985; Lamb et al. 1995; Ghosh & Lamb 1992). Stated differently, 
the magnetospheric beat-frequency model of the HBO requires that the neutron stars in the Z 
sources be near magnetic spin equilibrium. Indeed, the 200-350 Hz spin frequencies predicted 
by the model (see Ghosh & Lamb 1992) are much larger than the « 20-50 Hz HBO frequencies, 
as required. The similarity of the Z-source spin frequencies predicted by the magnetospheric 
beat-frequency model to the spin frequencies inferred from the separation frequencies of the 
kilohertz QPOs in the Z sources lends further support to the model (Miller et al. 1998) and to its 
implication that the neutron stars in the Z sources are near magnetic spin equilibrium (White & 
Zhang 1997; Psaltis & Lamb 1998). If they are, then their spin frequencies are given by (Ghosh & 
Lamb 1979, 1992) 



(Me)" 



(3) 
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where uJc is the critical fastness parameter and the angle brackets indicate an average over a time 
interval equal to the timescale on which the accretion torque changes the spin. 



Combining equations (|2[) and (^) and identifying i^hbo with z^mbf, we find 



t'HBO 



+ 1 



1 



(4) 



Equation (^) shows that the magnetospheric beat-frequency model of the HBO predicts that 
^HBo/^s + 1 should be ~ l/wc and hence only slightly larger than unity. The Z sources are thought 
to be accreting at near-Eddington accretion rates (see Lamb 1989; Hasinger & van der Klis 1989). 
The inner accretion disk in these sources is therefore expected to be radiation-pressure-dominated, 
in which case a « 0.2 (see Table 1). Hence, if depends only weakly on the instantaneous 
accretion rate M, then I'uBo/'^s + 1 will also depend only weakly on M and possibly also on the 
magnetic field and mass of the neutron star, through the dependence of lOc on these quantities 
(any such dependence is expected to be weak). 

According to the magnetospheric beat-frequency model, the HBO frequency i^bo is related 
to the frequency 1^2 of the upper kilohertz QPO only indirectly, through the dependence of both 
frequencies on the mass accretion rate. In all sources in which kilohertz QPOs have so far been 
discovered, 1^2 increases with inferred mass accretion rate (see, e.g., van der Klis et al. 1996; 
Strohmayer et al. 1996; van der Klis 1998). Here we explore the consequences of the simple 
ansatz 1^2 = i^o{M / Me)^ , where vq and A are constants that are specific to each source and may 
depend on the mass and magnetic field strength of the neutron star. This relation, with A ~ 1, is 
consistent with kilohertz QPO observations of several atoll sources, provided that the observed 
countrate from an atoll source is proportional to the mass accretion rate (see, e.g., Strohmayer et 
al. 1996; Ford et al. 1997). In all the Z sources, 1^2 is consistent with being ~ 1200 Hz when they 
are accreting at near-Eddington rates, which implies vq ~ 1200 Hz, independent of the expected 
modest differences in the masses and magnetic field strengths of these neutron stars. Using this 
simple ansatz, equation (0) can be written 



A a/A 



where 



and equation (^) becomes 



where 



^ + 1 = A,-H V-o)"/^ 

An = LOr. 



The inferred value of A2 therefore provides an estimate of the critical fastness uJc- 



(5) 

(6) 
(7) 

(8) 
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Fig. 1. — Correlations between vubo + i^s, the sum of the HBO frequency and the inferred spin frequency, 
and the upper kilohertz QPO frequency in five Z sources. Only a single typical error bar is shown for each 
source, for clarity. The lines show the best fits of the magnetospheric beat-frequency relation @ to the data. 
The error bars for the Sco X-1 points are smaller than the circles plotted. The error bars for the Cyg X-2 
points are substantially larger than for the other points, because of the large uncertainty in the Cyg X-2 
spin frequency. 

In order to test the relations (|5|) and (|^) predicted by the magnetospheric beat-frequency 
model, simultaneous measurements of i^hbO; z^2; and Vs are needed. The HBO and kilohertz 
QPO frequencies are directly observed, but oscillations at the neutron star spin frequency have 
not yet been detected in the persistent emission of any Z source. However, comparisons of the 
frequencies of the two simultaneous kilohertz QPOs observed in the persistent emission of the 
atoll sources with the frequencies of the nearly coherent oscillations observed in these sources 
during type I X-ray bursts indicate that the neutron star spin frequency is nearly equal to the 
frequency separation between the two kilohertz QPOs (Strohmayer et al. 1996; Miller et al. 
1998; Strohmayer et al. 1998; Psaltis et al. 1998; Miller 1999). Hence, for GX 17-^2, GX 5-1, 
GX 340-1-0, and Cyg X-2 we set the spin frequency equal to the average frequency separation of 
the kilohertz QPOs. In Sco X-1, the frequency separation of the kilohertz QPOs is consistent 
with being constant at the lowest inferred accretion rates but decreases at higher rates (van der 
Klis et al. 1997). Sco X-1 is thought to be accreting at near-Eddington mass accretion rates when 
the frequency separation of the kilohertz QPOs decreases, and it is therefore plausible that this 
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Fig. 2. — Correlation between i^hbo/^^s + 1 and the frequency of the upper kilohertz QPO in five Z sources. 
The sohd line is the best fit of the magnetospheric beat-frequency relation to the points with upper 
kilohertz QPO frequencies less than 850 Hz. 



decrease is related to the effects of radiation forces on the dynamics of the accretion flow near the 
neutron star (see, however, Mendez et al. 1998; Psaltis et al. 1998). Hence, in plotting the Sco X-1 
data, we set the spin frequency equal to the nearly constant frequency separation of its kilohertz 
QPOs at low inferred mass accretion rates. The spin frequencies we have adopted are listed in 
Table 2. 

Relation (|5|) describes adequately (x^ ^ 1-5 per degree of freedom) the dependence of the 
sum of the HBO and inferred spin frequencies on the frequency of the upper kilohertz QPO in the 
five Z sources in our sample, considered separately. Table 2 lists the best-fit parameters and their 
Icr errors and Figure 1 compares the best-fit relations with the frequency data on each source. 
If A 1, the power-law index a for all sources except Sco X-1 is !^ 0.2. This value is consistent 
with the expectation that the Z sources are accreting at near-Eddington rates and hence that 
the inner accretion disk is optically thick and radiation-pressure-dominated (see Table 1). For 
Sco X-1, however, a significantly weaker dependence of the HBO frequency on M is required, if A 
is independent of the mass accretion rate (but see Appendix A and Psaltis, Belloni, & van der 
Klis 1999). 

Figure 2 shows the quantity z^hbo/i^s + 1 plotted as a function of the upper kilohertz QPO 
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Fig. 3. — Confidence contours for tlie coefficient A2 and the power-law index a/A in the beat-frequency 
relation obtained by fitting this relation to the Z-source frequency data with upper kilohertz QPO 
frequencies less than 850 Hz and to the data of Sco X-1. The inner and outer contours show the 68% and 
99% confidence limits, respectively. The star indicates the best-fit values. The best-fit relation is compared 
with the frequency data in Fig. 2. Within the beat-frequency model of the HBO, the value of A2 is an upper 
bound on the critical fastness Ur- 



frequency 1^2 for the five Z sources in our sample. The frequency data on all the sources are 
consistent with a single, universal relation between z^hbOj ^s, and V2, as predicted by equation (|^, 
when V2 is < 850 Hz. This relation is shown as a solid line in Figure 2. Figure 3 shows the 
confidence contours for the power-law index a/ A and the coefficient A2 in relation (0) obtained 
by fitting this relation to all the data with V2 < 850 Hz. Assuming that i/q is ~ 1200 Hz, the 
best-fit value of A2 gives a lower bound on the critical fastness uJc, because is expected to be 
a monotonically increasing function of M and hence ((^M)") < (^Me)°. If the magnetospheric 
beat-frequency model is the correct explanation of the HBO, then uJc is ^ 0.8 for the magnetic 
field strengths and accretion rates of the Z sources. 

When f2 is > 850 Hz, the HBO frequencies of GX 17-1-2 are up to 2% higher than predicted 
by extrapolating the universal relation that holds at lower frequencies, whereas those of Sco X-1 
are as much as 5% lower. This indicates that there is at least one other important parameter that 
varies with i'2- For example, the structure of the inner disk may change at high accretion rates. 
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causing the exponent A to vary from source to source. This conjecture cannot be tested without a 
specific model for the variation in A, because if A is chosen to reproduce the behavior of the data, 
relation looses all predictive power. 

The magnetospheric beat-frequency model of the HBO requires that the neutron stars in 
all the Z sources be near magnetic spin equilibrium. The tight, universal correlation between 
the HBO, spin, and upper kilohertz QPO frequencies in all the Z sources when 1^2 < 850 Hz is 
explained by the model if A21'q^^ = uJcI^q^^ {{CM)") /{^Me)" is approximately the same in all of 
them. All the Z sources are thought to be accreting at very similar rates (comparable to the 
Eddington critical accretion rate; see, e.g.. Lamb 1989; Hasinger & van der Klis 1989) and hence 
((^M)")/(^Me)" is not expected to differ much from one to another. The critical fastness uJc 
is expected to be comparable to unity and may be a universal constant for a given inner disk 
structure (Ghosh & Lamb 1979, 1992). If so, then i^q^^ is nearly the same in the five Z sources in 
our sample and hence the frequency of the upper kilohertz QPO is a good absolute measure of the 
mass accretion rate in these sources. Stated differently, the relation between the upper kilohertz 
QPO frequency and the accretion rate appears to be very similar in all the Z sources. 

2.2. Discussion 

The analysis presented in §2.1 demonstrates that if (a) the neutron stars in the Z sources 
are spinning near their magnetic spin equilibrium rates, (b) the frequency separation between 
the upper and lower kilohertz QPOs is approximately equal to the neutron-star spin frequency, 
(c) the inner accretion disk in the Z sources is optically thick and radiation-pressure dominated, 
and (d) the upper kilohertz QPO frequency is proportional to the mass accretion rate through the 
inner disk, then the magnetospheric beat-frequency model is consistent with the observed behavior 
of the HBO frequency and, in particular, the tight, universal correlation (Figure 2) between the 
HBO and kilohertz QPO frequencies in all the Z sources in our sample. All of these assumptions 
are expected to be satisfied, as discussed in §2.1. 

As noted earlier, the magnetospheric beat-frequency model of the HBO is consistent with 
models for the upper kilohertz QPO that identify its frequency with an orbital frequency in the 
inner disk only if a small fraction of the accreting matter does not couple to the stellar magnetic 
field at the radius but instead remains in a geometrically thin Keplerian disk down to the radius 
responsible for the upper kilohertz QPO (i.e., ^ must be less than unity; see van der Klis 1998; 
Miller et al. 1998; and Alpar & Yilmaz 1997 for discussions). In particular, the interpretation of 
both the HBO and the two simultaneous kilohertz QPOs as rotational beat phenomena requires 
that there be two distinct radii in the inner accretion disk at which beating of the neutron star 
spin frequency with the orbital frequency produces a QPO, as is the case, for example, in the 
sonic-point model (Miller et al. 1998). 

Assuming that the HBO is a magnetospheric beat-frequency phenomenon, we can use the 



- 11 - 



observed HBO properties together with a general relation for the coupling radius to constrain 
the magnetic dipole moments of the neutron stars in the Z sources in a way that is largely 
independent of the structure of the inner accretion disk. In the Ghosh & Lamb (1979) model of 
disk-magnetosphere interaction, the radius at which the stellar magnetic field strongly couples 
to the gas in the accretion disk is given implicitly by (see Ghosh & Lamb 1991) 

4 \ 1/7 




for any model of the inner accretion disk. Here B^/Bp is the mean azimuthal magnetic pitch in 
the annulus of radial width Ar/r^ in the inner disk where the stellar field strongly interacts with 
gas in the disk. 

If the stellar magnetic field is too weak, it cannot couple strongly to the gas in the accretion 
flow well above the stellar surface and hence cannot generate magnetospheric beat-frequency 
oscillations. Hence in the magnetospheric beat-frequency model, the coupling radius must be 
larger than the neutron star radius i?NS) which requires 

where Mmax is the maximum mass accretion rate at which the HBO is detected. In deriving 



inequality (10) we have neglected the contributions of any higher multipole moments of the stellar 
magnetic field that may be present near the neutron star surface or may be induced by the 
electrical currents flowing in the disk (see Psaltis, Lamb, & Zylstra 1996 for a discussion). For the 
Keplerian frequency at the coupling radius to exceed the neutron star spin frequency, which is also 
required in the magnetospheric beat-frequency model (see Ghosh & Lamb 1979), the magnetic 
dipole moment must satisfy 



Bp I V 0.01 J V 300 Hz 



where we have used the fact that ^ Me. These upper and lower bounds on the magnetic 
dipole moment depend only very weakly on the neutron star mass. Figure 4 shows the resulting 
lower (eq. |jl^) and upper (eq. |11]) bounds on the magnetic dipole moments of the neutron stars 



in the Z sources, as a function of the relative width Ar/r^ of the coupling region. 

We can obtain an estimate of the magnetic dipole moment of the Z sources by using 
the value of Ai obtained by fitting relation (^) to the frequency data and the optically-thick, 
radiation-pressure-dominated model of the inner disk. The result is 

. (0.8-1.0) ( T^V' ( , (12) 



1200 Hz/ \2Mp, 
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Fig. 4. — Bounds on the magnetic moments of the neutron stars in the Z sources, derived from the 
magnetospheric beat-frequency model of the HBO. Solid lines: Upper bounds on the magnetic dipole 
moment derived from the requirement that the orbital frequency at the coupling radius be greater than 
the neutron star spin frequency. Dashed lines: Lower bounds on the magnetic dipole moment derived from 
the requirement that the coupling radius be larger than the radius of the star. Solid and dashed lines 
are plotted for two values (1 and 5) of the mean azimuthal magnetic pitch in the coupling region. Dotted 
line: Lower bound on the fractional radial width of the coupling region derived from the requirement that 
the predicted FWHM of the HBO be consistent with the observed width. Shaded region: Range of dipole 
magnetic moments estimated from the fits shown in Fig. 1 of the magnetospheric beat-frequency relation (^) 
to the Z-source data. The polar magnetic field strength that corresponds to a given magnetic moment is 
shown on the right vertical axis. For simplicity, the neutron stars in the Z sources were all assumed to have 
masses, radii, and spin frequencies of 2Mq, 10^ cm, and 300 Hz. 



where ^ may depend on the magnetic field strength. Relation (12) shows that the HBO frequencies 
predicted by the magnetospheric beat-frequency model are consistent with the HBO frequencies 
observed if /i27 ~ 1, which implies that the dipole magnetic fields of the Z sources have field 
strengths at the magnetic poles of ~ 10^ G for a 10 km neutron star. (Note that the estimated 
dipole magnetic moment depends only weakly on the unknown parameters ^ and vq.) 

The relative width Ar/r^ of the annulus where the stellar magnetic field strongly couples to 
the gas in the disk is expected to be greater than ~ 0.01 (Ghosh &; Lamb 1992). Its value can be 
bounded below using the observed FWHM of the HBO (see also Alpar & Shaham 1985; Lamb et 
al. 1985). Assuming that all other QPO broadening mechanisms — such as lifetime broadening — are 
negligible, we can estimate the relative width of the annulus in the accretion disk in which the 
interaction at the beat frequency affects the X-ray luminosity, from the relative width of the HBO 
peak in power spectra. The width of this annulus is necessarily smaller than the width Ar of the 
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layer where the magnetic field strongly interacts with the gas in the disk, and hence 

rm 3\UYlBO + l'sJ \10 Hz/ Vt'HBO + z^s/ 

where SiyuBO is the FWHM of the HBO. Figure 4 displays this bound on Ar/rm and the constraint 
it imposes on the dipole moment of the stellar magnetic field. 

Figure 4 shows that these additional physical bounds on the dipolar magnetic fields of the Z 
sources derived from the magnetospheric beat-frequency interpretation of the HBO are consistent 
both with each other and with the field strengths ~ 10^ G estimated in equation ( [T^ ) and by 
modeling the X-ray spectra of the Z sources (Psaltis et al. 1995; Psaltis & Lamb 1998). 



3. THE LENSE-THIRRING PRECESSION INTERPRETATION 

3.1. Model Predictions and Comparison with Observations 

In the nodal (Lense-Thirring and tidal) precession model of the HBO (Stella & Vietri 1998), 
a narrow ring or clumps of gas are assumed to be in a tilted orbit at the radius responsible for the 
upper kilohertz QPO and to precess with the frequency of a test particle in such an orbit (Stella 
& Vietri 1998). Alternatively, if the disk ends at this radius, one of the many weakly damped 
global precession modes localized near the inner edge of the disk (Markovic & Lamb 1998) may be 
excited. 

In the weak field limit, the nodal precession frequency of an infinitesimally tilted orbit at the 
radius where the frequency of a circular Keplerian orbit is i/k is (see Stella & Vietri 1998; Morsink 
& Stella 1999) 



(14) 



j/NP ~ 13.2 



M ) V300 Hzj U kHz/ 



Vm5/3; VO.Ol/ V300 Hzy VlkHz 



where rj = — (^//2^)(i/s/300 Hz)~^ in terms of A, the coefficient of the quadrupole moment of 
the gravitational field, and Izz, the neutron star moment of inertia with respect to its spin axis. 
Equation ( |l^ is derived by expanding the full expression for v^-p in a power series in Vg and 
retaining only terms up to second order. 

Lense-Thirring precession. — If the effects of the quadrupole component of the star's 
gravitational field are negligible, the localized warping modes of the inner disk will precess with a 
frequency close to the Lense-Thirring frequency of a test particle (Markovic &: Lamb 1998), which 
is given by the first term in equation (p!4[). Identifying the centroid frequency of the HBO with the 
Lense-Thirring frequency at the radius where the orbital frequency is equal to the frequency V2 of 
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Fig. 5. — Correlations between HBO frequency and upper kilohertz QPO frequency in five Z sources. The 
solid lines show the best fits of the Lense-Thirring relation (111) to the data. The dotted lines show the best 



fits of the more general power-law relation (Al), which is discussed in Appendix A. The HBO frequencies 
and lines for all the sources except GX 17+2 have been shifted downward by successive factors of two for 
clarity. 



the upper kilohertz QPO gives (Stella & Vietri 1998) 



mBO = -TT7l'sT^2 = 13.2 — 



M J V300 UzJ VI kHz 



1^2 



Hz 



(15) 



where / = 10*^^/45 g cm^ is the moment of inertia of the neutron star and 1^2 is the orbital 
frequency at the radius responsible for the upper kilohertz QPO. The X-ray visibility as well as 
the excitation and damping of the precession modes of the inner disk have not yet been addressed 
(see Markovic &: Lamb 1998 for a discussion). Equation (|l5|) predicts a relation between the HBO 
frequency, the spin frequency of the neutron star, and the frequency f2 of the upper kilohertz 
QPO that depends only on the structure of the neutron star, through the ratio I/M. 

Figure 5 shows the HBO frequencies observed in the five Z sources in the present sample, 
plotted against the frequencies of their upper kilohertz QPOs. Separate fits of equation ( p^ 
to the data on the individual sources, using the neutron star spin frequency inferred from the 
frequency separation of the kilohertz QPOs and treating I/M as a free parameter, give values of 
per degree of freedom of order unity for four of the Z sources but ~ 4 for Sco X-1. There is 
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Fig. 6. — Correlation between the HBO frequency i^hbo divided by the inferred neutron star spin frequency 
i^s (Table 2) and the upper kilohertz QPO frequency in five Z sources. Only a single typical error bar is 
shown for each source, for clarity. The solid line is the best fit of the Lense-Thirring relation (15) to all 
points with upper kilohertz QPO frequencies less than 850 Hz. The dotted line is the best fit of the more 
general power-law relation (^2) discussed in Appendix A to the same points. 



no other freedom in relation (15) and hence pure Lense-Thirring precession must be rejected as 
an explanation of the HBO. Moreover, as Stella & Vietri (1998) noticed (see also Wijnands et al. 
1998a; Jonker et al. 1998), the coefficients of i/^z^l required to fit the data give values of /45/M ^ 4, 
which is ~ 4-5 times larger than the largest ratios given by realistic equations of state for stars 
of any mass and about 2.5 times larger even than the largest ratio given by the extremely stiff 
relativistic mean-field equation of state L (see Table 3). 

In the Lense-Thirring precession model of the HBO, the relation between I'uBo/t^s and the 
upper kilohertz QPO frequency depends only on the mass of the star and the equation of state 



of neutron-star matter (see eq. |15|). Data from similar neutron stars should therefore follow 
similar relations. Figure 6 shows how the ratio I'uBO I'^s scales with the frequency vi of the upper 
kilohertz QPO; according to the Lense-Thirring precession model, this ratio should scale as v\. 
Figure 6 shows that the data with frequencies V2 < 850 Hz are consistent (Xdof ~ -'-•-'-) with a 



single relation of the form (15). Again, however, the coefficient given by the fit requires neutron 



stars with I45/M ^ 4, which is implausibly large. Furthermore, the points that have 1^2 > 850 Hz 
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Fig. 7. — Correlation in Sco X-1 between i^jbo/j^I ^^'^ ^^1^ ^ where z^hbo is the HBO frequency and V2 is 
the frequency of the upper kilohertz QPO. The solid line is the straight line that best fits the data. 

are inconsistent with the relation of the form ( |l5[ ) that fits the points with lower values of V2. 

Effect of classical precession. — Stella & Vietri (1998; see also the extended discussion in 
Morsink & Stella 1999) suggested that the flattening of the yuBO^i^2 correlation at high U2 might 
be caused by the increasing importance, as the disk penetrates closer to the star, of the classical 
precession caused the rotation-induced quadrupole component of the star's gravitational field. This 
precession is retrograde but smaller than the prograde gravitomagnetic precession and therefore 
tends to reduce the nodal precession frequency. We can test this suggestion quantitatively using 
the data for Sco X-1 and GX 17+2, which deviate most strongly from the relation of the form ( |l5|) 
that fits the points with low values V2 < 850 Hz.. 

If the HBO is caused by nodal precession and classical precession is important, then i^hbo/^^I 
should decrease linearly with increasing v^^^ (cf. Stella & Vietri 1998), because relation (|T^ ) can 
be rewritten as 

VI Hz/ VI kHz/ VMyV300 Hzy 

"^■^ (ii^) (om) (soolk) (n^) ■ ^^^^ 
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Figure 7 plots against 1^2 for Sco X-1 and also shows the best-fit straight line with slope 

1/3, which has Xd o f ~ ^- '^'^^ ^^^^ ^^^^ data in Figure 7 can be fit satisfactorily by a straight 
line with slope 1/3 is not strong evidence for this scaling, because the range of measured 1^2^^ 
values is very narrow. However, we can use the best-fit value of the intercept of the straight line 
with the vertical axis to estimate the value of the parameter rj that characterizes the quadrupole 
moment of the gravitational field and to estimate I/M. The results are 

,Seo = 2^3xlO-M=/3(^)" (17) 

and /r \ / \ 

/J45W^^\ (18) 
\M J V300 Hz/ ^ ^ 

The value oi I/M required by equation (|l^) is ~ 20 time larger than the largest values given 
by realistic neutron-star equations of state (see Table 3). The deviation of the GX 17+2 data 
from the power-law relation (|15|) at high frequencies (see Fig. 6) requires that the classical 
precession frequency be negligible for 1^2 < 850 Hz but comparable to the Lense-Thirring precession 
frequency at slightly larger values of z^2- This is not possible, because the exponents of 1^2 in 
the Lense-Thirring and classical precession terms of equation (jlj) are too similar. The data are 
therefore inconsistent with the predictions of the simple nodal precession model. 



3.2. Discussion 

In the Lense-Thirring precession model of the HBO, the HBO frequency, the spin frequency 



of the neutron star, and the frequency of the upper kilohertz QPO are related by equation (15). 
As discussed in the previous section, the data are consistent with this relation when the frequency 
of the upper kilohertz QPO is < 850 Hz, but the inferred value of I/M is implausibly large (see 
also Stella & Vietri 1997, 1998; Morsink & Stella 1999). Aside from the rather unhkely possibility 
that neutron stars have ratios of I/M that are four times as large as the largest values for stellar 
models constructed with realistic equations of state, three other possibilities have been suggested 
for reducing this large discrepancy. 

First, the frequency difference Au between the kilohertz QPOs might be equal to half the 
neutron star spin frequency Vs rather than equal to it. This is very unlikely in any beat-frequency 
model of the kilohertz QPOs, because it would require a special direction that rotates with the 
neutron star but affects the inner accretion disk only once every two beat periods. However, 
Av = i^fs appeared possible, given the initial analysis of the data taken during the type I X-ray 
bursts of 4U 1636—536 (Zhang et al. 1997; Strohmayer et al. 1998), which showed a strong 
oscillation at about 580 Hz, approximately twice the frequency separation of the two kilohertz 
QPOs. However, further analysis of this data by Miller (1999) using a matched-waveform filtering 
technique has revealed the presence of a weak coherent oscillation at about 290 Hz, approximately 
equal to the frequency separation of the two kilohertz QPOs. Thus, it now appears very unlikely 
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that the spin frequencies of these neutron stars are twice the frequency separation of their kilohertz 
QPOs. 

Second, the observed HBO frequencies and their second harmonics might represent the 
second and fourth harmonics of the fundamental Lense-Thirring frequency (|l5|), rather than the 
first and second harmonics. Indeed, a precessing circular orbit has a two- fold symmetry that 
could, in principle, produce even-order harmonics that are stronger than the odd-order harmonics. 
Moreover, power-density spectra of the Z sources show a relatively strong, broad-band noise 
component at frequencies comparable to the ones predicted by relation ([l5|). This so-called 
low-frequency noise (Hasinger & van der Klis 1989) might inhibit detection of the fundamental 
of a low- frequency precession frequency (see, e.g.. Fig. 6a of Kuulkers et al. [1994] for peaked 
features in the low-frequency noise component of GX 5—1). Determination of the upper limits 
on the amplitudes of any QPOs at these frequencies would significantly constrain this possibility. 
Note, however, that even if the HBO and its overtone are the second and fourth harmonics of the 
precession frequency, the I /M ratios required to explain the HBO observations would still be a 
factor ;^ 2 larger than predicted by realistic neutron-star equations of state. 

A further difficulty with the Lense-Thirring precession interpretation of the HBO is that the 
observed correlation between the HBO and kilohertz QPO frequencies is significantly different from 
what is predicted by relation (|l5| ) when the frequency of the upper kilohertz QPO is > 850 Hz. 
As demonstrated in §3.1, this difference cannot be explained by classical precession; nor can it be 



explained by strong- field corrections to relation (15) (Stella & Vietri 1998) 



A third possibility is that radiation forces increase the ratio of z^np to 1^2 by the factor ~ 2-5 
required to bring it into agreement with the observed HBO and upper kHz QPO frequencies. 
The Z sources are thought to be accreting at near-Eddington mass accretion rates when the 
frequencies of the upper kilohertz QPOs are comparable to ~ 1 kHz (see, e.g., Psaltis et al. 1995). 
Hence radiation forces, which were neglected in equation (|l4|), almost certainly are important. At 
near-critical luminosities, both orbital and nodal precession frequencies can be altered by large 
factors compared to their values in the absence of radiation; hence radiation forces might possibly 
explain the large discrepancy between the observed frequencies of the HBO and the frequencies 
predicted by the nodal precession model. 

An explanation in terms of the combined effects of Lense-Thirring precession and radiation 
forces would, however, require the physically implausible result that radiation forces leave the 
variation with V2 basically unchanged while increasing the ratio of t'NP to 1^2 by a factor ~ 2-5. 
Such an explanation would also require that the QPO peaks not be significantly broadened by the 
radiation drag force at the same time that radiation forces are strong enough to change the orbital 
and precession frequencies by a factor ~ 2-5. The HBO peaks in Sco X-1, for example, have 
fractional widths Sv/iy <^ 0.5 even when the inferred accretion rate is near the Eddington critical 
rate (van der Klis et al. 1997). More fundamentally, if radiation forces do change the orbital and 
precession frequencies of gas accreting onto the Z sources by large factors, as they may well do, the 
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observed correlation between the HBO and upper kilohertz QPO frequencies would be explained 
primarily by the effect of the radiation forces and not by the gravitomagnetic torque. 

4. CONCLUSIONS 

In §2 and §3 we have studied in detail the behavior of the HBO frequencies observed in five Z 
sources and in particular their correlation with the frequencies of the kilohertz QPOs, comparing 
the observed behavior with the behaviors predicted by the magnetospheric beat-frequency (Alpar 
& Shaham 1985; Lamb et al. 1985) and Lense-Thirring precession (Stella k, Vietri 1998) models of 
the HBO. 

In §2 we showed that the magnetospheric beat-frequency model is consistent with the 
observed correlation between the HBO and upper kilohertz QPO frequencies in the five Z sources 
studied here if, as expected, the neutron stars in these sources are spinning near their magnetic 
spin equilibrium rates, the frequency separation between the upper and lower kilohertz QPOs is 
approximately equal to the neutron-star spin frequency, the inner part of their accretion disks are 
optically thick and radiation-pressure-dominated, and the frequency of the upper kilohertz QPO 
is approximately proportional to the mass accretion rate. The model predicts a universal relation 
between the horizontal branch oscillation, stellar spin, and upper kilohertz QPO frequencies that 
agrees well with the data on five Z sources. The spin rates predicted by the model are consistent 
with the range of the spin frequencies of the Z sources inferred from the frequency separation of 
their kilohertz QPOs if they are all accreting at similar, near-critical rates and all have 10^-10^^ G 
dipole magnetic fields. Such magnetic fields are consistent with models of Z-source X-ray spectra. 
The inferred value of the critical fastness for the accretion rates and magnetic field strengths of 
the Z sources is ^ 0.8. If the frequency of the upper kilohertz QPO is an orbital frequency in the 
accretion disk, the magnetospheric beat-frequency model requires that a fraction of the accreting 
gas does not couple strongly to the stellar magnetic field until it has penetrated to within a few 
kilometers of the neutron star surface. 

In §3, wc showed that the trend of the correlation between the HBO frequency and the upper 
kilohertz QPO frequency observed at upper kilohertz QPO frequencies V2 < 850 Hz agrees with 
the trend predicted by the Lense-Thirring precession model. However, the observed trend is 
inconsistent with the model for U2 > 850 Hz. The observed magnitudes of the HBO frequencies 
are ^ 4-5 times smaller than the magnitudes predicted by the Lense-Thirring precession model 
for realistic neutron-star equations of state. Thus, in order to be consistent with the observed 
magnitudes, either I/M must be ^ 4-5 times larger than expected or the principal frequency of 
the X-ray oscillation generated by nodal precession must be ^ 4-5 times the nodal precession 
frequency. 

We thank Greg Cook for making available the numerical code used to compute the suite 
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APPENDIX 



A. AN EMPIRICAL DESCRIPTION OF THE CORRELATION OBSERVED 
BETW^EEN THE HBO AND UPPER KILOHERTZ QPO FREQUENCIES 

In this appendix we show that the correlations observed between the HBO frequency z^hbO; 
the upper kilohcrtz QPO frequency f2, and the frequency separation Ai^ between the two kilohertz 
QPOs can be characterized by simple power-law relations among the frequencies involved. Our 
purpose is to facilitate comparison of the present data with future data (see, e.g., Psaltis et al. 
1999) or other theoretical models. 

The frequency correlations in the five Z sources in our sample can be described adequately by 
the power-law relation 

vmio = 13.2 ai (y^) ' ' 

where the constants ai and hi are different for each source. The confidence contours obtained by 
the fitting this relation to the measured HBO and upper kilohertz QPO frequencies of the five 
Z sources in the present sample are shown in Figure 8. The power-law index that describes the 
Sco X-1 data is significantly smaller than the index that describes the data on the other four 
sources. The best-fit relations for each source are the dashed lines shown in Figure 5. 

For all the Z sources except Sco X-1, the best-fit value of the parameter oi is approximately 

proportional to the spin frequency inferred from the frequency separation of the two kilohertz 
QPOs (sec Fig. 6). Indeed, when the upper kilohertz QPO frequency is < 850 Hz, the correlation 
between i^jbOj and 1^2 is described adequately by the relation 

'^~ = i'-^»^fe)(r^)^ (^^' 

with 02 ~ 4.6, and 62 ~ 1-8. The confidence contours obtained by fitting this relation to the 
Sco X-1 points and to the points on the other four sources for which 1^2 < 850 Hz are shown in 
Figure 9. As Figure 6 shows, the frequency correlation is significantly flatter when 1^2 is > 850 Hz. 
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Confidence contours for the coefficient a\ and the power-law index 6i, obtained by fitting 



Fig. S 

equation (Al) to the frequency data on each of five Z sources, considered individuaUy. The inner and 
outer contours show the 68% and 99% confidence limits, respectively, while the stars indicate the best-fit 
values for each source. The best-fit relations are shown as dashed lines in Fig. 5. 



In Sco X-1, the HBO and kilohertz QPOs were simultaneously detected mostly when it was on 
the normal branch. In the other four Z sources, the HBO and kilohertz QPOs were simultaneously 
detected mostly when they were on their horizontal branches. The transition from the horizontal 
to the normal branch is thought to take place when the mass accretion rate increases to within 
a few percent of the Eddington critical rate (see Lamb 1989; Psaltis et al. 1995). If so, the 
resulting change in the accretion flow pattern (see Lamb 1989) might be responsible for the 
different dependences of the HBO frequency z^hbo and the instantaneous frequency separation Ai/ 
of the two kilohertz QPOs on the upper kilohertz QPO frequency V2 in Sco X-1, compared to the 
dependences in the other four Z sources in our sample. 

The ratio of i^ibo to Ai/ in Sco X-1 increases more steeply with than does the ratio of fneo 
to the (constant) inferred spin frequency Vg. This is demonstrated most clearly by the correlation 
plots shown in Figure 10. (For all the Z sources in the sample except Sco X-1, Ai^ is consistent 
with being constant or with varying in the same way as it does in Sco X-1 [Psaltis et al. 1998; see 
also Wijnands et al. 1997, 1998b; Jonker et al. 1998]. This is true mostly because of the relatively 
large uncertainties in the measured kilohertz QPO frequencies.) The dependence of vy^bqI IS.v on 
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Fig. 9. — Confidence contours for tlie coefficients 02 and power-law indices 62 derived by fitting relation ( A2) 
to two disjoint sets of frequency data. The contours marked "Sco X-1" were obtained by fitting this relation 
to the Sco X-1 data. The contours marked 'V2 < 850 Hz" were obtained by fitting this relation to the 
frequency data on all the other sources, when the upper kilohertz QPO frequency is less than 850 Hz; the 
best-fit version of this relation is shown by the dashed line in Fig. 6. The inner and outer contours show the 
68% and 99% confidence limits, respectively, while the stars indicate the best-fit values for each source. 



1/2 in Sco X-1 is more consistent with the behavior seen in the other four Z sources and suggests 
that we plot i/hbo/^^^ against 1^2 for all five of the Z sources in our sample. The result is shown in 
Figure 11. In this plot we have included only points derived from simultaneous observations of the 
HBO and kilohertz QPO frequencies. The larger scatter of the points in Figure 11 at 1^2 < 850 Hz 
compared to the scatter of the points in Figure 7 corresponding to the same values of i^hbo and 
1^2 is caused by the large uncertainties in Av. The points plotted in Figure 11 are consistent with 
power-law relations of the form 

b3 



fHBO = 13.2 03 



1^2 



(A3) 



.300 RzJ VI kHz, 

The confidence contours obtained by fitting this relation to the data on each of the sources except 
Cyg X-2 are shown in Figure 12. Cyg X-2 was not included because its HBO and kilohertz QPO 
frequencies were measured simultaneously only once. These contours show that a single universal 
relation of the form (A3) with 03 = 4.2 and 63 = 1.6 is consistent with the data on all the sources 
except Sco X-1. In fact, 63 = 1.6 is consistent with all the data. The Sco X-1 and GX 17-1-2 data 
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Fig. 10. — Top panel: Correlation 
observed in Sco X-1 between j^hbo/'^s, 
the HBO frequency divided by the 
neutron star spin 

frequency inferred from the frequency 
separation of the two kilohertz QPOs 
at low countrates, and the frequency 
of the upper kilohertz QPO. Bottom 
panel: Correlation between i^HBo/^i^i 
the HBO frequency divided by the 
instantaneous frequency separation of 
the kilohertz QPOs, and the frequency 
of the upper kilohertz QPO. In both 
panels the dashed line shows the best 



fit of relation (A.3) to the data 
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have relatively small uncertainties and give best-fit coefficients 03 that are slightly but significantly 
different. The uncertainties in the data on GX 5—1 and GX 340+0 are sufficiently large that their 
contours allow values of 03 that are consistent with either the Sco X-1 or the GX 17+2 value. 

The surprising correlation between the HBO frequency, the frequency of the upper kilohertz 
QPO, and the instantaneous frequency separation of the kilohertz QPOs shown in Figure 11 
may be coincidental. The relatively large uncertainties in the currently measured kilohertz QPO 
frequencies in all the Z sources except Sco X-1 prevent us from drawing any firm conclusions 
about the significance of this correlation. However, the strikingly similar relation between the 
lower and upper kilohertz QPO frequencies in all LMXBs that show kilohertz QPOs (Psaltis et 
al. 1998) together with the correlation shown in Figure 11 suggests that the varying frequency 
separation of the kilohertz QPOs in Sco X-1 is a general property of the kilohertz QPOs and is 
related, directly or indirectly, to the frequency of the HBO. Additional data are needed to test 
directly this conjecture (see Psaltis et al. 1999 for an alternative possibility). 
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Fig. 11. — Correlation between the HBO frequency j^hbo divided by the instantaneous frequency separation 
Av of the kilohertz QPOs and the frequency of the upper kilohertz QPO in five Z sources. The uncertainties 
in the GX 17+2 and Sco X-1 points are typically smaller than the sizes of the symbols plotted for these 
sources and have therefore been omitted. The dashed line is the best fit of relation (^^) to the frequency 
data on all the sources except Sco X-1. 
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Table 1. Scaling of the Keplerian Frequency at the Coupling Radius^ 



Disk Model'' 


VKfi (Hz) 


a 


/3 


7 


IG 


430 


0.38 


-0.87 


0.82 


IR 


210 


0.23 


-0.77 


0.70 


2B 


80 


0.72 


-0.86 


0.43 


2S 


50 


2.55 


-1.20 


-0.60 



'^Prom Gliosh & Lamb 1992. 

'^IG: Optically tliick, gas-pressure dominated (GPD) disk; IR: Optically thick, 

radiation-pressure dominated (RPD) disk; 2B: Two-tcmpcraturc. optically thin 
GPD disk with Comptonized bremsstrahlung; 2S: Two-tcmpcraturc, optically 
thin GPD disk with Comptonized soft photons (for references to these disk models 
see Ghosh & Lamb 1992). 



Table 2. Best-Fit Parameters for the Magnetospheric Beat-Frequency 

Model 



Source 


i^s (Hz)** 


Normalization (^i) 


Exponent (a/ A) 


GX 17+2 


293.5±5.5 


365±1 


0.21±0.01 


GX 5-1 


312.5±9.2 


378±3 


0.16±0.02 


GX 340-hO 


324.5±9.3 


389±4 


0.15±0.03 


Cyg X-2 


324.4±44.9 


393±1 


0.17±0.02 


Sco X-1 


304.7±1.3 


354±1 


0.06±0.03 



'^Inferred from the frequency separation of the kilohertz QPOs (see text). 
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Table 3. Maximum I/M Ratios for Illustrative Neutron Star Equations 

of State^ 





Mass (Mq) 


/45/M 


A 


1.66 


0.63 


AU 


2.13 


0.93 


FPS 


1.80 


0.76 


L 


2.70 


1.73 


M 


1.80 


0.91 



''The tabulated values are the maximum ratios found in a computational survey 
of all stable masses and spin rates for the equation of state indicated, using the 
code developed by Cook, Shapiro, & Teukolsky (1994). 

^A, AU, and FPS are realistic equations of state; M illustrates the effects of 
a strongly repulsive tensor interaction; L is an extremely stiff equation of state 
produced by a simple relativistic mean-field theory. For a recent discussion of 
these equations of state, see Miller, Lamb, & Cook (1998). 



